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Abstract 

Various quasi-exact solvability conditions, involving the parameters of the periodic associated 
Lame potential, are shown to emerge naturally in the quantum Hamilton-Jacobi approach. It 
is found that, the intrinsic nonlinearity of the Riccati type quantum Hamilton-Jacobi equation 
is primarily responsible for the surprisingly large number of allowed solvability conditions in the 
associated Lame case. We also study the singularity structure of the quantum momentum function, 
which yields the band edge eigenvalues and eigenfunctions. 
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I. INTRODUCTION 



Quasi-exactly solvable (QES) Hamiltonians, interconnecting a diverse array of physical 
problems, have been the subject of extensive study in recent times jll2l3l4l5IU] . These 
systems, containing a finite number of exactly obtainable eigenstates, have been linked with 
classical electrostatic problems, as also to the finite dimensional irreducible representations of 
certain algebras. Some of these studies employ group theoretical methods, others are based 
on the symmetry of the relevant differential equations |2l3l4j . The key to the existence of the 
finite number of identifiable states is the quasi-exact solvability condition, relating certain 
potential parameters of these dynamical systems. An interesting feature, distinguishing these 
QES systems from known exactly solvable cases, is the presence of complex zeros of the wave 
functions. Hence, quantum Hamilton- Jacobi (QHJ) formalism, being naturally formulated 
in the complex domain, is ideally suited for studying the QES problems Although 
polynomial potentials have been studied rather exhaustively, QES periodic potentials have 
not received significant attention in the literature. 

The associated Lame potential (ALP) 

CIl ^ ( OC Tfi^ 

V{x) = a{a + l)msn^(x,m) + 6(6 + l)m — ^ — (1) 

dn (x,m) 

is an interesting example of a periodic potential, which is exactly solvable, when a = b and 
shows QES property, when a ^ b. Here, sn (x, m), cn (x, m) and dn(x, m) are the doubly 
periodic elliptic functions with modulus parameter m [7j. The ALP has a periodic lattice of 
period K{m) with the basis composed of two different atoms which are alternately placed. It 
possesses a surprisingly large variety of QES solvability conditions depending on the nature 
of the potential parameters a and 6. 

For a, 6 being unequal integers, with a > 6 > 0, there are a bound bands followed by a 
continuum band. If a — 6 is odd (even) integer, it has 6 doubly degenerate band edges of 
period 2K{4:K), which can not be obtained analytically. For a, 6 having half- integral values 
(with a > b), there are infinite number of bands with band edge wave functions having 
period 2K {4:K), if a — 6 is odd (even). Of these infinite number of bands, a — 6 bands have 
band edges, which are non-degenerate with period 2K{4:K) and 6 + | doubly degenerate 
states of period 2K{4K) which can be obtained analytically. For a being an integer and 6 
being a half integer or vice versa, one can obtain some exact analytical results for mid-band 
states IHlElIin]. 
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It is quite natural to enquire about the origin of this rich QES structure in the associated 
Lame potential |llll2j . This paper is devoted to the study of the same through the quantum 
Hamilton- Jacobi approach. The quasi-exact solvability conditions, involving the parameters 
of the periodic associated Lame potential, are shown to emerge naturally. The intrinsic 
nonlinearity of the Riccati type QHJ equation is responsible for these allowed solvability 
conditions in the associated Lame case. We also study the singularity structure of the 
quantum momentum function, which yields the band edge eigenvalues and eigenfunctions. 

In our earlier studies, we had looked at non-periodic ES, QES and ES periodic potentials 
through the QHJ formalism which was initiated by Leacock and Padgett fH^ IT^ . We 
were successful in obtaining the quasi exact solvability condition [5] for QES models and 
in obtaining the eigenvalues and eigenfunctions for both the ES ^H] and QES models jS]. 
Within the same formalism, we have also proposed a new method to solve ES periodic 
potentials [1^] and obtain the band edge eigenvalues and eigenfunctions. The present study 
makes use of the QHJ formalism as devolved in our earlier papers ^3 IHj and pH] • We refer 
the interested readers to them for all the details. 

The QHJ formalism revolves round the logarithmic derivative of the wave function p, 



known as the quantum momentum function (QMF), which satisfies a Riccati type equation 



We have found that, the knowledge of the singularity structure of the QMF and the residue 
at these singular points are the only information needed to obtain the QES condition and the 
solutions. In all our earlier studies, it was assumed that, the point at infinity is an isolated 
singularity, which turned out to be true in all the cases. With this same assumptions on p, 
we proceed to obtain the QES condition and the forms of the band edge wave functions for 
the general ALP in the next section. In this paper we concentrate on the cases where, a and 
b are either integers or half integers. In section HI, we analyze the situation where both a, 
b are integers taking, the values 2 and 1, respectively. The case when a and b are both half 
integers, with values 7/2 and 1/2 respectively, is analyzed in section IV. 



p = —ih— \nip 
ax 



(2) 



p- 




2M{E-V{x)). 



(3) 
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II. QES CONDITION AND THE FORMS OF THE WAVE FUNCTIONS 



The QHJ equation for the ALP, putting K = 2M = 1, is given by 

(cn ^ fx) \ 
E-a(a + l)msn^(x) -6(6+ l)m — oM" • (4) 
dn [X) J 

Note that with the transformation b —b— 1 or a — > —a — 1, the potential does not change. 
Hence, for our analysis, without loss of generality, we take a, b positive, with a > b. Defining 
p = —iq and substituting it in (jD), one obtains 

do cn ^ f 

+ — = a(a + l)m sn ^(x) + b(b + l)m — - ^- (5) 
ax dn (x) 



Changing the variable to t = sn (x), and writing 



q=^{l-t'){l-mt^)<P, with = x+-(^_-_^ + -_^j, (6) 

Equation (0), gets transformed into 

2 dx mH^ + 2m{l - 2b{b + 1)) 2 + t^ 2E - mt^{l - 2a{a + I)) _ 
^ ^~dt^ 4(l-mt2)2 ^ 4(1-^2)2 + 2{l-t^){l-mt^) ~ °' ^ ^ 

For all our later calculations, we shall treat x as the QMF and ((Tj) as the QHJ equation. 

Singularity Structure : From ((Tj), we see that x has fixed poles at t = ±1 and 
t = In addition to the fixed poles, the QMF has finite number of moving poles and no 

other singular points in the complex plane. Hence, one can write x a sum of the singular 
and analytical parts as follows 

bi b\ di di PL ^, ^ /„\ 

where bi,b[,di and d'l are the residues at t = ±1 and respectively, which need to be 

calculated. P„ is an n^'^ degree polynomial with -p^ = J2k=i being the summation of 
terms coming from the n moving poles with residue one. The function Q{t) is analytic and 
bounded at infinity. Hence, from Liouville's theorem it is a constant, say C. The residues 
at the fixed poles can be calculated by taking the Laurent expansion around each individual 
pole and substituting them in 0. Comparing the coefficients of different powers of t, one 
gets two values of residues at each pole owing to the quadratic nature of the QHJ equation. 
Thus the two values of the residues, 6i ,b[ at t = ±1 are 



At t = ±-^, one has 

3 6 16 ^.3616 ,^ . 

(ii = - + -, and rf' = - + -, . (10) 

^4 2' 42 ^4 2' 42 ^ ^ 

Since, there is no way of ruhng out one of the two values of these residues, we need to 
consider both the values. We demand h\ = b[ and di = d[, a condition coming from the 
parity constraint = — 

Behaviour at infinity : Equation (jH)) gives the behaviour of x in the entire complex 
plane. Hence, for large t, 

26i + 2rfi + n 

X ^ , (11) 

where the restriction 6i = b[ and di = d\ has been applied. This should match with the 
leading behaviour of x obtainable from ((Jj). Note that, the assumption on the singularity 
structure on x is equivalent to the point at infinity being an isolated singularity. Since x 
has at most an isolated singular point at infinity, one can expand x iii Laurent series around 
the point at infinity as, 

x(t) = Ao + ^ + ^ + --- . (12) 

Substituting (fT^ in ((Tj) and comparing various powers of t, one gets Aq = 0, thus making 
Q{t) in (jHl) equal to zero. Further one obtains 

Ai = a + 1, -a. (13) 

Since both the equations (|TT|l and (fT3|l give the leading behavior of x at infinity, both should 
be equal. Thus 

26i + 2rfi + n = Ai. (14) 

Taking various combinations of hi and di from Q and ()10j). substituting them in ()14|) one 
obtains the QES condition for each combination as given in table 1 for A = a + 1. Thus one 
sees that all the allowed combinations of residues give one of the forms of QES condition |H] , 
where n = 0, 1, 2 ■ ■ ■ . Note that, the other value of A, i.e., —a, when substituted instead of 
a+1 in ()14j) . gives the QES condition for negative values of a, 6 i. e., for a —>■ — a— 1, 6 —>■ —6—1 
in 6i, di. 

Forms of wave function : From one can write in terms of p as 

■ijj{x) = exp ipdx^ . (15) 



Changing the variable to t and writing p in terms of x? one gets 

*W^e.p(/(x4(^ + ^))*). (16) 

Substituting x from (jHJ in the above equation gives the wave function in terms of the residue 
61, di and the polynomial P„ : 

In terms of the original variable x the wave function takes the form, 

i){x) = (cnx)"(dnx)^P„(snx) (18) 

where a = P = Hence, for each set of 61, di one gets a wave function given by 

|T8|l . The degree n of this polynomial, which is obtained from (|Tl|l as, 

n = a + 1 - 26i - 2c/i (19) 

is in terms of either a + b or a — b, as evident from table 1. The forms of the wave function 
can be found and are given in table 1, for the two different cases, when a + b and a — b are 
odd and even separately. 

Case 1, Both a + b, a — b are even: We introduce = ^ and M = where M and 
are integers, and obtain the forms of the wave functions in table 2, in terms of M and A^, 
for the four sets of combinations of 61 and di in table 1. 

Case 2 Both a + b, a - b odd: Introducing A^' = ^ and M' = where M' and 
A^' are integers, we obtain the wave functions, in terms of M' and A^', for the four sets of 
combinations of bi and di in table 1. 

From the forms of the wave functions in tables 2 and 3, one observes that the number of 
linearly independent solutions is different for the two cases. The unknown polynomial in the 
wave function can be obtained by substituting x from (jH)) in the QHJ equation (jZj), which 
gives 

Pnit) + 4P„(t) + ^1^) + G(t)P„(t) = (20) 
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where 



_ t^Abj - 2h + \)-2b, + l mH\Ad\ - 2rfi + \) - 2md^ + m( ^-y 
~ if - 1)2 + (mt2 - 1)2 + 

2E + {IQhidi - 1 - 2a(a + l))mt 
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2(1 -t2)(l _mt2) 

The above differential equation is equivalent to a system of n linear equations for the coef- 
ficients of the different powers of t in Pn{t). The energy eigenvalues are obtained by setting 
the corresponding determinant equal to zero. In the next section we obtain the band edge 
wave functions for the associated Lame potential, when a = 2 and 6 = 1. 



III. ALP WITH a, b INTEGERS : 

For this case, we consider the associated Lame potential with a = 2, b = 1. For the 
purpose of comparison with literature, we work with the supersymmetric potential 

2 

CXI tJj 

y_(x) = 6msn (x) + 2m 2 " 4m (21) 

dnx 

This potential is same as with a = 2 and 6=1, except that a constant has been added 

to make the lowest energy equal to zero. The QHJ equation in terms of x is 

2 , m^t^ _ 6m 2 + f 2E + 8m - ISmt^ _ 
^ ^ ^ ^ 4(1 - mt2)2 + 4(1 - t2)2 + 2(1 - t2)(i _ ^^2) - 0- i ) 

Apart from n moving poles, x has poles at t = ±1 and t = -^Xj^Jm. As in the previous 
section one can write Xi with the parity constraint as 

This gives the form of x the entire complex plane, where P„ is yet to be determined. Note 
that for this potential the combination a + 6 and a — h are both odd ^.e, 3 and 1 respectively. 
Hence, we use table 3 to obtain all the information regarding the residues at the fixed poles, 
number of moving poles of Xi number of linearly independent solutions and their form, for 
each set, by taking the values of a = 2, 6 = 1, M' = and A^' = 1. The unknown polynomial 
in the wave function can be obtained from ()20|). where G{t) for this potential satisfies, 

_ t\Abl - 2b, + i) - 26i + I mH^jAdj - 2di + i) - 2mrfi - ^ 
~ (t2 - 1)2 + (rnt^ - 1)2 + 

2E + 8m + {IQbidi - 13)mt^ 



2(l-t2)(i_rnt2) 



(24) 



Using (pUj) and (j21I), one gets the exphcit expressions for the eigenfunctions and the eigen- 
values as given in table 4. From the table, we see that the first set of residues gives n = — 1, 
which will not be considered as n cannot be negative. Thus this particular case of Lame 
potential has five band edge solutions, which can be obtained analytically out of an infinite 
number of possible states. 



IV. ALP WITH a, b HALF INTEGERS : 

The potential studied here is the supersymmetric associated Lame potential, with a = 
7/2, 6=1/2: 

63 2 3 cn^x ^ 29 , 

V- = — msn X + -m — 2 ^ r"^ + "9 v^^J 

4 4 dn X 4 

where Sg = ^/4 — 4m + 25m2. In terms of x(t)y the QHJ equation takes the form, 

2 , 2 + mH'^ -m 4E + 8 + 29m - 4^9 - 65mt^ 

^ + ^ + 4(t2 - 1)2 + 4(mt2 _ ly + 4(t2-l)(mt2-l) ~ ' ^ ' 

Note that, for this case, a + b = 4 and a — 6 = 3, these are even and odd respectively. Hence, 
for such cases, one needs to use sets 1 and 3 from table 2 and sets 2 and 4 from table 3 
in order to get the four groups of the eigenfunctions. The solutions for this potential are 
given in table 5. We see that there is a degeneracy in the band edge energy eigenvalue, 
14 — 7m + 6g. All the solutions agree with the known results [H]. 



V. CONCLUSIONS 



In conclusion, in this study, we have demonstrated the applicability of QHJ formalism to 
QES periodic potentials. We have been successful in obtaining the quasi-exact solvability 
conditions and band edge solutions for cases when both a and b are integers or half-integers. 
The origin of the large number of solvability conditions for the ALP case, comes out naturally 
in QHJ approach. Interestingly, the singularity structure of the QMF for the QES periodic 
potentials is similar to that of the ES periodic potentials. This structure is completely 
different from those of the polynomial potentials. The case where one of the parameter a 
and b is an integer, the other being a half-integer, requires further careful study and will be 
reported elsewhere. 
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TABLE I: The quasi exact solvability condition from the four permitted combinations of 61 and di 
for the general associated Lame potential. 



set 


bi 


di 


2bi + 2di + n = Xi 


QES condition 


1 


3/4 


3,6 
4 2 


2+b+n= a 


b-a = -n -2 


2 


3/4 


1 b 
4 2 


1 — b + n = a 


a-l-b-l-1 = n-l-2 


3 


1/4 


3,6 
4 2 


1 + 6 + ra = a 


b-a = -n-1 


4 


1/4 


1 b 
4 2 


— 6 + n = a 


a+b = n 



TABLE II: The form of the wave functions for the four sets of residue combinations when a + b 
and a — b are even and equal to 2A'^ and 2M respectively. 



set 


bi 


di 


n = Ai - 26i - 2di 


n{M, N) 


wave function 'ip{x) 


LI solutions 


1 

2 
3 
4 


3/4 
3/4 
1/4 
1/4 


3,6 
4 2 

1 6 
4 2 

3,6 
4 + 2 

1 6 
4 2 


a-b-2 
a + b-1 

a-b-1 
a + b 


2M-2 
2N-1 

2M - 1 
2N 


cn x{dn ■^+^P2m-2 (sn x) 

(r.>^2iv-i(snx) 
idnx)^+^P2M-i{snx) 

P2N(snx) 
(dnx)'' 


M 
N 
M 
N + 1 



TABLE III: The form of the wave functions for the four sets of residue combinations when a + b 
and a — b are odd and equal to 2N' + 1 and 2M' + 1 respectively. 



set 


bi 


di 


n = Ai - 26i - 2di 


n{M',N') 


wave function ip{x) 


LI solutions 


1 
2 
3 
4 


3/4 
3/4 
1/4 
1/4 


3,6 
4 + 2 

1 6 
4 2 

3,6 
4 + 2 

1 6 
4 2 


a-6-2 
a + 6- 1 
a-6-1 
a + b 


2M' - 1 

2N' 
2M' 
2N' + 1 


cn x(dn a;)^+''P2M'-i (sn a:;) 

(-,P2A^,(snx) 

(dnx)''+iP2M'(snx) 

'P2Ar'+i(sna;) 
(dna;)'' 


M' 

N' + 1 
M' + 1 
N' + l 
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TABLE IV: For the Associated Lame potential V-{x) = 6m sn (x)^ + 2m^^-^|yj — 4m, with a = 2 



and 6=1 the residues, the value of n, number of linear independent solutions, the band edge 
eigenf unctions and eigenvalues are as follows. Here a + b = 3 and a — 6 = 1 which give N' = 1 and 
M' = 0. 



set 


bi 


di 


n 


LI solutions 


eigenfunction ■^(a;) 


eigenvalues 


1 


3/4 


5/4 


-1 








2 


3/4 


-1/4 


2 


2 


^ (3msn - 2 ± V4 - 3m) 


5 - 3m ± 2V4 - 3m 


3 


1/4 


5/4 





1 


dn^x 





4 


1/4 


-1/4 


2 


2 


^ (3msn ^x — 2 — m±\/4 — 5m + m2 


5 - 2m ± 2Vm'^ - 5m + 4 



TABLE V: For the Associated Lame potential F_(x) = ^msn^x + f - 2 - ^m + (Jg with 

a = 7/2 and 6=1/2 the residues, the value of n, number of linear independent solutions, the band 
edge eigenfunctions and eigenvalues are as follows. Here a + 6 = 4 and a — 6 = 3 which give N = 2 
and M' = 1. 



set 


bi 


di 


n 


LI solutions 


eigenfunction '^(x) 


eigenvalues 


1 


3/4 


1 


1 


1 


cnx(dnx)^/2snx 


^9 — m + 2 


2 


3/4 





3 


2 


cnx(dnx)^/2snx 


6g —m + 2 












cnx(dnx)~"'^/2snx(l — 2sn'^x) 


14 - 7m + Sg 


3 


1/4 


1 


2 


2 


(dnx)^/2(12msn2x - 5m - 2 - Jg) 















(dnx)^/2(12msn2x - 5m - 2 + 5g) 


2dg 


4 


1/4 





4 


3 


(dnx)^/2(12msn2x — 5m — 2 — 6q) 















(dnx)^/2(12msn2x - 5m - 2 + 5g) 


2Sg 












1 — Ssn^xcn^x 


14 - 7m + Sg 
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